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Abstract 

The process e + e~ — > 77(7) is studied using data recorded with the OPAL detector 
at LEP. The data sample taken at a centre-of-mass energy of 189 GeV corresponds to a 
total integrated luminosity of 178 pb _1 . The measured cross-section agrees well with the 
expectation from QED. A fit to the angular distribution is used to obtain improved limits 
at 95% CL on the QED cut-off parameters: A + > 304 GeV and A_ > 295 GeV as well as a 
mass limit for an excited electron, M e * > 306 GeV assuming equal e*e7 and ee7 couplings. 
Graviton exchange in the context of theories with higher dimensions is excluded for scales 
G+ < 660 GeV and G_ < 634 GeV. No evidence for resonance production is found in the 
invariant mass spectrum of photon pairs. Limits are obtained for the cross-section times 
branching ratio for a resonance decaying into two photons and produced in association 
with another photon. 
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1 Introduction 



This paper reports a study of the process e + e~ — > 77(7) using data recorded with the OPAL 
detector at LEP at an average centre-of-mass energy of 188.63 ± 0.04 GeV with an integrated 
luminosity of 178.3 ± 0.4 pb _1 . This is one of the few processes dominated by QED even at high 
LEP energies. Since the differential cross-section is well known from QED |], 0, any deviation 
from this expectation hints at non-standard physics processes contributing to these photonic 
final states. Any non-QED effects within the general framework of effective Lagrangian theory 
are expected to increase with centre-of-mass energy ||. A comparison of the measured photon 
angular distribution with the QED expectation can be used to place limits on the QED cut-off 
parameters A±, contact interactions (e + e~77) and non-standard e + e~7-couplings as described 
in Section 3. The possible existence of an excited electron e* which would change the angular 
distribution M, is investigated as well as the possibility of graviton exchange in the context 
of theories with higher dimensions 0. In addition, a search is made for the production of a 
resonance X via e + e~ — > X7 followed by the decay X — > 77, using the invariant mass spectrum 
of photon pairs in three-photon final states. The process e + e~ — > 77(7) has been analysed 
previously at lower energies f| [7], [|, |, |Hj. The selection presented in this paper is almost 
unchanged with respect to the previous OPAL analysis at y/s = 183 GeV 0. 



2 The OPAL detector and Monte Carlo samples 



A detailed description of the OPAL detector can be found in Ref. |TlJ. The polar angle 9 is 
measured with respect to the electron beam direction and is the azimuthal angle. For this 
analysis the most important detector component is the electromagnetic calorimeter (ECAL) 
which is divided into two parts, the barrel and the endcaps. The barrel covers polar angles 
with I cos#| < 0.82 and consists of 9440 lead-glass blocks in a quasi-pointing geometry. The 
two endcaps cover the polar angle range 0.81 < | cos^l < 0.98 and each consists of 1132 blocks. 
For beam-energy photons, the spatial resolution is about 11 mm, corresponding to 0.2° in 9, 
and the energy resolution is about 2% in the barrel and 3 - 5% in the endcaps, depending on 
the polar angle. The ECAL surrounds the tracking chambers. Hit information from the central 
jet chamber CJ and the central vertex drift chamber CV is used to reject events which are 
consistent with having charged particles coming from the interaction point. Outside the ECAL 
is the hadronic calorimeter HCAL which is incorporated into the magnet yoke, and beyond that 
are the muon chambers. Both the HCAL and the muon chambers are used to reject cosmic ray 
events. 

Various Monte Carlo samples are used to study the selection efficiency and expected background 
contributions. For the signal process e + e~ — > 77(7) the RADCOR |2j generator is used while 
FGAM [|12|] is used for e + e~ — > 7777. FGAM does not take into account the electron mass 
leading to a divergent cross-section at small angles and hence can not be used if at least one 
photon is along the beam axis. The Bhabha process is simulated using BHWIDE and 



TEEGG [fLl ]. The process e + e — > 1/1/7(7) * s simulated using KORALZ JT^]. All samples were 



processed through the OPAL detector simulation program [HJ and reconstructed in the same 
way as the data. 
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3 Cross-section for the process e + e 



77 



The Born-level differential cross-section for the process e + e 



77 in the relativistic limit of 



lowest order QED is given by |T7 



da 



) 



Born 



Oi 



y 2 1 + cos 2 
s 1 — cos 2 



(1) 



where s denotes the square of the centre-of-mass energy, a is the fine-structure constant and 
is the polar angle of one photon. Since the two photons are identical particles, the event angle 
is defined by convention such that cos is positive. 

In a multi-photon event, it is important to choose an appropriate definition of the event angle. 
The event angle is not uniquely defined since the two highest-energy photons in general are not 
exactly back-to-back. The event angle cos 0* used in this paper is defined as 



where 0\ and 02 are the polar angles of the two highest-energy photons. This definition was 
chosen such that the deviations in the angular distribution with respect to the Born level are 
small, in this case between 3 and 6% for cos#* < 0.99, as was shown in Ref. 0. For two-photon 
final states, cos0* is identical to | cos#| and for three-photon events in which the third photon is 
along the beam direction, 0* is equivalent to the scattering angle in the centre-of-mass system 
of the two observed photons. Since the angular definition is based on the two highest-energy 
photons, events with one of these escaping detection along the beam axis are rejected from 
the analysis. The observed angular distribution is corrected for these higher-order effects using 
0(a 3 ) Monte Carlo studies in order that it can be directly compared to the following model 
expectations given in lowest order. 

In Ref. 0, possible deviations from the QED cross-section for Bhabha and M0ller scattering 
are parametrised in terms of cut-off parameters A±. These parameters correspond to a short- 
range exponential term added to the Coulomb potential. This ansatz leads to modifications of 
the photon angular distribution of the form 



Alternatively, in terms of effective Lagrangian theory, a gauge-invariant operator may be added 
to QED. Depending on the dimension of the operator, different deviations from QED can be 
formulated [|J. Contact interactions (77e + e~) or non-standard 7e + e~ couplings described by 
operators of dimension 6, 7 or 8 lead to angular distributions with different mass scales A. 
In most cases these deviations are functionally similar ||. Therefore, only the cross-section 
predicted by a dimension-7 Lagrangian, given by 



is studied here. Limits on mass scales for Lagrangians of dimension 6 and 8 can easily be 
derived from limits on A± and A' ||. 






(3) 




(4) 
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The existence of an excited electron e* with an e*e7 coupling would contribute to the photon 
production process via t-channel exchange. The resulting deviation from ( ^ ) depends on 

V""/Born 

the e* mass M e » and the coupling constant k of the e*e7 vertex relative to the ee7 vertex ||: 

= f ) + /(M e *,K,s,cos0). (5) 




dtt J 



Born 



The function / is explicitly given in Ref. 0. In the limit M e » ^> y/s, (Jjj^) t approaches (Jj^ 
with the mass related to the cut-off parameter by M e * = s/k A+. 



A± 



Recent theories || have pointed out that the graviton G might propagate in a higher-dimensional 
space where additional dimensions are compactified while other Standard Model particles are 
confined to the usual 3+1 space-time dimensions. This would allow the large Planck scale Mpi 
with an effective scale M D of order of the electroweak scale ((9(10 2_3 GeV)) via Mpj = R n M£ +2 , 
where R is the compactification radius for n higher dimensions. Although the contribution from 
e + e~ — > G* — > 77 for a single mode is small compared to the Standard Model contribution, 
the very large number of possible excitation modes in the extra dimensions might lead to a 
measurable effect on the resulting differential cross-section, given by 





± — G ± 4 (1 + cos 2 6) + -f- G ± s (1 - cos 4 9) , (6) 



7T 



64tt 2 



G± Born 

where G± are scale parameters of order of the effective scale M-q. 



4 Event selection 



A) Preselection 

The preselection has only small changes compared to our previous analysis H. Events are 
selected if they have two or more ECAL clusters within the polar-angle range | cos#| < 0.97 
each with more than 1 GeV of deposited energy uncorrected for possible energy loss in the 
material before the ECAL. Other energies used in the analysis are corrected for these losses. A 
cluster must consist of at least two lead-glass blocks. This preselection is tightened by requiring 
a total energy in the ECAL of at least 10% of the centre-of-mass energy. 



B) Rejection of non-physics backgrounds 

A cosmic-ray particle can pass through the hadronic and electromagnetic calorimeters without 
producing a reconstructed track in the central tracking chambers. The requirements to reject 
these events are the same as in our previous analysis. We reject events with 3 or more track 
segments found in the muon chambers. In the case of fewer than three such track segments, the 
event rejection depends on the highest-energy HCAL cluster of the event. Events are rejected if 
this HCAL cluster is separated from each of the photon candidates by more than 10° in azimuth 
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and has at least 1 GeV of deposited energy in the case of one or two muon-track segments, or 
at least 15% of the observed ECAL energy if no muon-track segments are found. 

Another type of background is consistent with well localised electronic noise in the ECAL. 
The cuts used to reject these events are slightly relaxed compared to our previous analysis to 
account for the larger cluster size due to the higher centre-of-mass energy. An event is rejected 
if one of these accumulations consists of more than 14 ECAL clusters or has an extent of more 
than 0.5 rad in 9 or 0.5 rad/sin# in <fi. 

C) Kinematic requirements 

The selection is based primarily on the requirement of small missing energy and small missing 
transverse momentum and is unchanged with respect to our previous analysis. Non-physics 
background and events containing invisible particles or having only a small energy deposit in 
the ECAL are rejected by this selection. In particular, events from the Standard Model process 
e + e~ — > 1/1/77 fail this selection. Better resolution on quantities such as the total longitudinal or 
transverse momentum of an event is obtained by using cluster angles and three-body kinematics 
where possible, rather than using measured energies. Therefore the event sample is divided into 
four classes, I- IV. The classes are distinguished by the number of photon candidates iV 7 , the 
acollinearity angle £ aco i, and the aplanarity angle £ a pia n : 

£acol = 180° -ai2 (7) 

faplan = 360° - (a 12 + (X 13 + CX 23 ) , (8) 

where is the angle between clusters i and j and the indices are ordered by decreasing cluster 
energy. 

All events having £ acol < 10° (i.e. with the two highest-energy clusters almost collinear) are 
assigned to class /, independent of the number of photon candidates. More than 90% of the 
events fall in this class. Acollinear events (£ aco i > 10°) are separated into three classes. Events 
with N y — 2 are assigned to class II, typically containing an energetic photon that escapes 



Event class 


I 


II 


III 


N 


Number of 
photon candidates 

Acollinearity 
Aplanarity 


Ny > 2 

£acol < 10° 


AT 7 = 2 

£acol > 10° 


Ny = 3 

£acol > 10° 
Kaplan < 0.1° 


Ny > 3 

£acol > 10° 
Kaplan > 0.1° 


Energy sum 
Transverse momentum 
Longitudinal momentum 


El > 0.6^/i 


Eg > 0.6v/i 

B < 0.2 
-Eiost < E 1 ,E 2 


Ef > 0.6^ 
Jh < O.ly/S 

Pi < E 1 ,E 2 


Ef > 0.6^ 

p t < 0.1 v/i 
Pi < E 1 , E 2 



Table 1: Summary of the kinematic cuts. The upper part of the table describes the definition of 
the four event classes, the lower part describes the applied cuts. For definition of the variables 
see the text. In the case of more than three observed photons in class IV there is no requirement 
on the aplanarity. 
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detection near the beam axis (|cos#| > 0.97). Planar events (£ a pian < 0.1°) with exactly 3 
observed photon candidates are assigned to class III. The remaining events, i.e. aplanar three- 
photon events (£ ap ian > 0.1°) and events with more than three observed photon candidates are 
assigned to class IV. Classes III and i\^are distinguished for two reasons. Firstly, the main 
signal Monte Carlo generator calculates only events with up to three photons, class IV events 
are not described by this program. Secondly, for class III events, three-body kinematics can 
be used to calculate the photon energies from the angles, to determine the invariant masses of 
photon pairs. This is not possible for class IV events. The event class definitions are summarised 
in Table |I|. 

Events in all classes are required to meet the following criteria: 

CI The energy sum of all observed ECAL clusters must be more than 60% of the centre-of- 
mass energy and the transverse momentum of the event must be less than 10% of the 
centre-of-mass energy. 

C2 The longitudinal momentum of the event must be smaller than the energy of each of the 
two highest-energy clusters to ensure that the two highest-energy photons in the event 
are observed. 

For class I events it is assumed that missing energy is negligible due to the collinearity. There 
is only one cut applied to the sum of the two highest cluster energies E s = E x + E 2 . The 
distribution of Eg/y/s is shown in Figure |I[ Events are selected if Eg > 0.6y/s. 

For class 77 events three-body kinematics can be used to calculate missing energy assuming 
a third unobserved photon with |cos#| = 1. The missing longitudinal momentum can be 
estimated as E\ ost = a/s/ [1 + (sin^i + sin^)/] sin (6*i + 6* 2 )|] and the transverse momentum 
can be approximated by the imbalance B = (sin#i +sin02) |cos [(0i — 2 )/2]|, where 6\, 6 2 and 
0i, 02 are the angles of the observed clusters. The energy sum Eg = E\ + E 2 + Ey ost is calculated 
by summing the two observed cluster energies and the missing energy. Events are selected if 
Eg > 0.6a/s, B < 0.2 and E\ ost is less than both Ei and E 2 . This last requirement ensures 
that the two highest-energy photons are those observed. Figure 0a shows the distribution of 
Eg / yfs for the data before and after the cut on B, compared with the signal Monte Carlo. 

For events in classes III and IV, the cluster energies must be used in addition to the cluster 
angles in calculating the transverse and longitudinal momenta (p t , p\) of the system. The 
energy sum E s u = Y^2\ Ei + p\ is calculated by adding p\ to the cluster energies Events are 
selected if Eg > 0.6\/s and p t < 0.1y/s. Again, the missing energy along the beam axis, now 
determined by p\, must be smaller than the energies of the two highest-energy clusters. Figure 
shows the distribution of Eg 1 J y/s for the data before and after the cut on p t , compared with 
the signal Monte Carlo. 

The non-physics background is reduced to a negligible level after the above kinematic require- 
ments which are summarised in Table |l|. 
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D) Charged event rejection 

Bhabha events, for example, have electromagnetic cluster characteristics similar to 77(7) events, 
but can be distinguished by the presence of tracks in the central tracking chambers. The 
rejection of all events having tracks in the central tracking chambers, however, would lead to 
an efficiency loss because of photon conversions. Nevertheless, contributions from any channel 
with primary charged tracks must be reduced to a negligible level for this analysis. To achieve 
this, for each cluster, hits in the inner part of the drift chambers CV and CJ which are associated 
in <p with the cluster are counted. If the number of such hits in CV (CJ) is above a ^-dependent 
threshold, as described in 0, a CV (CJ) token is assigned to the cluster. There are three ways 
to reject an event: 

1. The single veto rejects events with both a CV and a CJ token assigned to the same cluster. 

2. The double veto rejects events with either a CV or a CJ token assigned to the two 
highest-energy clusters. 

3. Events are rejected if there is a reconstructed track separated by more than 10° in azimuth 
from all photon candidates. 



5 Corrections and systematic errors 

Table ^] shows the number of selected events from the data, the signal Monte Carlo and the 
major background sources after the different selection cuts. Other processes contribute only 
a small number of events [|7|. The preselected data have some contribution from non-physics 
backgrounds and from the process e + e~ — > qq until the kinematic cuts are applied. There is 
little efficiency loss up to this point. The restriction on the missing longitudinal momentum 
(-Ejost for class il) rejects events with a high-energy photon escaping along the beam axis. 
About 97% of the remaining sample consists of Bhabha events and is well described by the 



Cut 


Data MC 77(7) e + e (7) e7(e) z/z/77 


A) Preselection 

B) Non-physics bg. reject. 
CI) Kinematic requirements 
C2) E loBt , Pl <E 1 ,E 2 

D) Charged event rejection 


116 129 90 576 2 297 82 157 6 078 44 
96 307 90 393 2 296 82 014 6 039 44 
88 715 87 429 2 273 79 848 5 305 3.2 
78 782 79 460 2 109 77 183 168 0.05 
1 740 1 777 1 775 0.4 0.9 0.04 


Statistical errors 


42 8 8 0.2 0.4 0.01 



Table 2: Number of selected events after the different cuts described in the text. The numbers 
are given for the data and the sum of Monte Carlo samples with the breakdown by final states 
given in the following columns. The BHWIDE generator is labelled by e + e~(7) and TEEGG 
by e7(e). All Monte Carlo predictions are normalised to the integrated luminosity of the data. 
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Monte Carlo. The principal decrease in efficiency is due to the charged event rejection. After 
the final selection 1.3±0.5 background events are expected which is negligible compared to the 
expected 77(7) signal. 

The efficiency and angular resolution of the reconstruction are determined using an O (a 3 ) 
Monte Carlo sample with full detector simulation. In the angular range cos#* < 0.8 the 
efficiency is in the range 91 - 95% and drops to 40 % for cos#* > 0.96. The requirement with 
the lowest efficiency is the charged event rejection leading to a loss of events with early photon 
conversions. 

To compare the conversion probability in data and Monte Carlo, a special event sample is 
selected. This sample consists of events with exactly 2 clusters, an acollinearity angle £ aco i < 5° 
and an energy sum E I s >Qfiy/s. At least one cluster must be consistent with a non-converted 
photon: neither a CV or a CJ token, nor a track with a minimum number of hits must be 
assigned to the cluster. This leads to a sample of 2810 unbiased clusters, where the cluster is 
selected via the opposite photon, i.e. all clusters, where the opposite cluster is consistent with 
a non-converted photon. This sample is used to study the probability that a CV or CJ token 
is assigned to a cluster. There are 13.5 ± 1.5 background events expected from Bhabha and 
multi-hadronic events in the angular range cos 6* < 0.8. This may be compared to 14 events 
found by scanning the data after removing an additional 8 events originating from radiative 
return to J/^. Concerning the single-veto probability, 51 clusters with both a CV and a CJ 
token assigned are found in the data, in agreement with the 77(7) Monte Carlo expectation of 
44 clusters. Also for the double-veto probability for events that are not rejected by the single 
veto, a good agreement between data and Monte Carlo is observed: 41 clusters with either a 
CV or a CJ token assigned are found in the data while 40 are expected from Monte Carlo in 
the angular range cos#* <0.8. In the angular range cos8* <0.97 there are 258 clusters observed 
in the data and 248 expected from Monte Carlo. The study of the data alone leads to a 0.5% 
systematic error per cluster on the single- veto probability at cos 9 = due to the data statistics 
only. Including the information from the Monte Carlo, the overall systematic error on the 
efficiency is taken to be 1% at cos# = 0, corresponding to 20% of the inefficiency; this value is 
used for all angles. The 1% error is taken to be correlated between all cos^* bins. 

The veto probability due to detector occupancy is determined from randomly triggered events. 
It is found to be 0.7% for the charged-event rejection (D) and 0.3% for the non-physics back- 
ground rejection (B). 

The agreement between generated and reconstructed event angles gives an event-angle reso- 
lution of about 0.2°. In addition, the cluster angle has been compared to the track angle for 
Bhabha events. For clusters with | cos#| > 0.96, the cluster angle is systematically about 0.4° 
closer to the beam axis than the track angle. For clusters with | cos#| < 0.94, the difference is 
less than 0.1°. Due to the cut-off at cos9* < 0.97 this would lead to a decrease of the measured 
total cross-section by 1.1% if this effect is caused by the cluster angle. It is included in the 
systematic error on the total cross-section. 

The luminosity is derived from small-angle Bhabha scattering measured on both sides of the 
detector in the polar angle region 34 mrad < 9 < 56 mrad. Uncertainties derive from the 
selection and the theoretical cross-section, as well as from a 20 MeV uncertainty on the beam 
energy. To be conservative and include additional uncertainties an increased systematic error 
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of 0.5% is taken into account for the final errors. 

Since the deviations from QED (Equations § - ^) are given with respect to the Born level, the 
observed angular distribution is corrected to the Born level. The effect of radiative corrections 
is quantified by the ratio 1Z of the angular distribution of the e + e~ — ► 77(7) Monte Carlo and 
the Born-level cross-section: 

It is assumed that the higher-order corrections in the context of the studied models are equal 
to those expected from QED. The ratio TZ varies between 1.03 and 1.06 within the studied 
angular range cos#* < 0.97 and is used to correct the data bin by bin to the Born level. A 1% 
systematic error from 0{a A ) and higher-order effects is taken to be correlated between all bins. 



6 Results 

Cross-section 

In Table [3] the numbers of observed events in the different kinematic classes are compared to the 
Monte Carlo expectation after the final selection. Since there is no Monte Carlo generator for 
the case of four-photon events where at least one photon is along the beam axis, no prediction 
for class TV events with three observed photons is given. The prediction for events with four or 
more observed photons is calculated using FGAM. One class TV event with 4 detected photons 
is observed. Figure [I] shows a class / event with 4 detected photons, where the lowest-energy 
photon has an energy of 13.5 GeV. The expectation for a class / event with four photons with 
at least 10 GeV each is 0.07 events. The total cross-section cr for the process e + e~ — > 77(7), 
determined from 1740 events selected in the range cos#* < 0.97, is also given in Table [3| The 
cross-section is corrected for detection efficiency, 0(a 3 ) effects derived from TZ (Eq. ||) and the 
veto probability due to detector occupancy. The cross-section at the different LEP energies as 
measured by OPAL in the range cos#* < 0.9 is shown in Figure |j. All measurements are in 
good agreement with the QED prediction. 



Class 


/ 


II 


III 


TV 


All 


Cross-section [pb] 




> 2 


2 


3 


3 


> 4 


events 


cos#* < 0.9 


cosfl* < 0.97 


Observed 


1587 


125 


20 


7 


1 


1740 


7.55 ±0.18 ±0.14 


11.4 ±0.3 ±0.3 


Expected 


1612 


131 


32 




0.8 


1776 


7.49 


11.8 



Table 3: Number of observed and expected events in the angular range cos 6* < 0.97 for different 
classes. AT 7 is the number of observed photon candidates. No Monte Carlo expectation is 
available for class TV events with 3 observed photons. The total cross-section is corrected to 
the Born level and is given for two angular ranges to allow comparison of these results with 
those of previous OPAL measurements. The first error is statistical, the second is systematic. 
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The angular distribution of the observed events and the measured differential cross-section 
obtained by applying efficiency and radiative corrections, are shown in Figure |5|. The data have 
a x 2 / NDF = 14.8/20 with respect to the QED expectation (solid line). The 95% CL interval of 
a x 2- fit f° the data of the function (|§) A (Eq. |3|) is also shown (dashed line). To obtain the 

limits at 95% confidence level, the probability is normalised to the physically allowed region, 
i.e. A + > and A_ < as described in Ref. For both functions (|§) A and (^) A , the 

X 2 distribution is parabolic as a function of the chosen fitting parameters A^ 4 and A' -6 . The 
X 2 distribution for (|§) G from graviton exchange in the context of extra dimensions is also 

approximately parabolic. The asymmetric limits xf 5 on the fitting parameter are obtained by: 



So 9S T(x,fi,a)dx 
/ °° T(x, /i, a)dx 



0.95 



and 



I°-r(x,/i,o-)dx 



0.95 



(10) 



where T is a Gaussian with the central value and error of the fit result denoted by fi and 
o", respectively. This is equivalent to the integration of a Gaussian probability function as 
a function of the fit parameter. It is not possible to integrate the probability function as a 
function of A±, A' or G±. Therefore the limits on the fit parameters are derived via Eq. |i~0| . 
The 95 % CL limits on the cut-off parameters are derived from the limits on the fit parameters, 
e.g. the limit on G + is obtained as [^(G^ 4 )] -1 / 4 . 

The limit on the mass of an excited electron M e * as a function of the coupling-constant ratio k 
for the (e*e7)-vertex, which is fixed during the fit, is shown in Figure ^. In the case of (jgi) 

the cross-section does not depend linearly on the chosen fitting parameter M~ 2 and the limit 
corresponds to an increase of the x 2 by 3.84 with respect to the minimum. 

The fit results are summarised in Table |j. Limits obtained from a combined fit to the data 
taken at 189 and 183 GeV are also given. The limits obtained are 30-115 GeV higher than 
our previous results M. 



Fit parameter 


Fit result 


(189 GeV) 


95% CL Limit [GeV] 
189 GeV 183+189 GeV 






10- 12 GeV" 4 


A+ 


> 


297 


304 




(-36 ±71) ■ 












A_ 


> 


281 


295 


A'- 6 


(-2.7 ±6.3) 


• 10~ 18 GeV" 6 


A' 


> 


674 


672 




(-"IS) ■ 


10~ 12 GeV" 4 


G + 


> 


660 


660 


Gt 


G- 












> 


600 


634 


M- 2 


(0.0 ±6.6) ■ 


10- 6 GeV" 2 




> 


311 


306 



Table 4: Fit results and 95% CL lower limits obtained from the fit to the differential cross- 
section. The last column shows the limit obtained from a combined fit to the data taken at 
183 and 189 GeV. The limit for the mass of an excited electron is determined assuming the 
coupling-constant ratio k — 1. 
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Resonance production 



A resonance X produced by the process e + e~ — > X7 and decaying into two photons, X — > 77, 
would be visible in the two-photon invariant mass spectrum, since this process leads to a three- 
photon final state without missing energy. Searches for such a resonance have been performed 
previously at the Z° peak and at higher energies 0, ||, ^] , leading to bounds on Higgs and 
gauge boson interactions |£J. For this search, 20 events from class III are used. The invariant 
mass of each photon pair is shown in Figure ^a). The energies of the three photons are not 
based on the measured cluster energies but are calculated from the photon angles assuming 
three-body kinematics: 

E k oc sin aiij ; E\ + E 2 + E 3 = y/s, (11) 

with E k the energy of one photon and ctij the angle between the other two photons. In this 
case, the typical mass resolution for photon pairs is about 0.5 GeV. Deriving the mass from the 
measured cluster energies would lead to a resolution of about 3%. The distribution is consistent 
with the Monte Carlo expectation from the QED process e + e~ — > 77(7). There is no evidence 
for an enhancement due to a resonance. 



An upper limit on the total production cross-section multiplied by the two-photon branching 
ratio is calculated using the method of Ref. p4] . This method uses fractional event counting 
where the weights assigned to each photon pair depend on the expected resolution and the 
difference between the hypothetical and the reconstructed mass. The limits shown in Figure [7]b) 
are obtained assuming the natural width of the resonance to be negligible. The e + e~ — > 77(7) 
background is subtracted. For the efficiency correction, the production and subsequent decay 
of the resonance are assumed to be isotropic. The mass range is limited by the acollinearity 
restriction. Regarding a model with anomalous couplings of the Higgs boson p3] , this analysis 
gives access to a larger mass range than the analysis of e + e~ — > HZ with H — > 77 [25 . However, 
in the region of overlap the direct search is more sensitive. 



7 Conclusions 



The process e + e~ — > 77(7) has been studied using data taken with the OPAL detector at centre- 
of-mass energies up to 189 GeV. The measured angular distribution and total cross-section for 
this process both agree well with QED predictions. The limits (95% CL) on cut-off parameters 
are A + > 304 GeV, A_ > 295 GeV and A' > 672 GeV. An excited electron is excluded for 
M e * < 306 GeV assuming the e*e7 and ee7 coupling to be the same. Graviton exchange in 
the context of theories including extra dimensions is excluded for scales G + < 660 GeV and 
G- < 634 GeV. In the 77 invariant mass spectrum of events with three final-state photons, no 
evidence is found for a resonance X decaying to 77. A limit on the production cross-section 
times branching ratio is derived as a function of the mass Mx- 
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OPAL 




Figure 1: Scaled energy sum E^/y/s for class / events, a) shows events with cos#* < 0.8 and 
b) shows events with cos6>* > 0.8. The points show the data and the histogram represents 
the Monte Carlo expectation, normalised to the integrated luminosity of the data. The cut is 
indicated by an arrow. 
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Figure 2: a) Scaled energy sum E 1 * / \/s for class II events after the cut on Ei ost . b) Scaled 
energy sum E 1 / 1 / ' y/s for class III and W events after the cut on p\. The solid points show the 
data after the cut on the imbalance B (Fig. a) and the transverse momentum p t (Fig. b), and 
the open points show the data suppressed by these requirements. The histogram represents 
the Monte Carlo expectation, normalised to the integrated luminosity of the data. The Monte 
Carlo distribution is shown after the cuts with have negligible effect. The cut is indicated by 
an arrow. 
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Side v i ew - plane of Thrust axis 




E1 =82.7GeV 
E2 =45.2 GeV 
E3 -29.5 GeV 
E4=13.5 GeV 



ECAL 



HCAL 
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Figure 3: Display of an event with four detected photons. Reconstructed clusters in the ECAL 
are shown as blocks where the size is proportional to the energy deposit. Since the two highest 
energy photons have an acollinearity of 3.9° it is selected as class I event. There are 0.07 events 
expected in class / with four photons with an energy of more than 10 GeV each. 
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Figure 4: Total cross-section for the process e + e~ — > 77(7) with cos#* < 0.9. The curve corre- 
sponds to the Born-level QED expectation. The lower plot shows the distribution normalised 
to the QED expectation. The data are corrected for efficiency loss and higher-order effects 
and correspond to the Born level. Results at lower energies are taken from previous OPAL 
publications || 0, |l9t1 . 
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Figure 5: The measured angular distribution for the process e + e~ — > 77(7) at ^/s = 189 GeV. 
The histogram shows the observed number of events per bin. Note the smaller width of the 
highest cos 9* bin. The points show the number of events corrected for efficiency and radiative 
effects. The inner error-bars correspond to the statistical error and the outer error-bars to the 
total error. The solid curve corresponds to the Born-level QED prediction. The dashed lines 
represent the 95% CL interval of the fit to the function (j^) 
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Figure 6: Lower limit (95% CL) on the mass M e * of an excited electron as a function of the 
square of the e*e7 coupling constant ratio k 2 obtained from a combined fit to data taken at 
183 and 189 GeV. 
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Figure 7: Results of a search for resonance production in class III events, a) shows the invariant 
mass of photon pairs for data (points) and for the e + e~ — > 77(7) Monte Carlo expectation (open 
histogram). The hatched histogram represents a 77 resonance at 95.5 GeV with a cross-section 
times branching ratio of 0.05 pb. The binning is chosen to match the expected mass resolution, 
b) shows the upper limit (95% CL) for the cross-section times branching ratio for the process 
e + e~ — > X7, X — > 77 as a function of the mass of the resonance X. 
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